Evidence for nodal superconductivity in LaFePO 
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In several iron-arsenide superconductors there is strong evidence for a fully gapped supercon- 
ducting state consistent with either a conventional s-wave symmetry or an unusual s± state where 
there the gap changes sign between the electron and hole Fermi surface sheets. Here we report 
measurements of the penetration depth A(r) in very clean samples of the related iron- phosphide 
superconductor, LaFePO, at temperatures down to ~ 100 mK. We find that A(r) varies almost 
perfectly linearly with T strongly suggesting the presence of gap nodes in this compound. Taken 
together with other data, this suggests the gap function may not be generic to all pnictide super- 
conductors. 



A key question in understanding the physics of the 
recently discovered iron pnictide superconductors is the 
origin of the pairing interaction. Measurements of the 
symmetry and anisotropy of the superconducting energy 
gap are very useful for helping to decide between com- 
peting theories. Although several experiment designed to 
deduce the gap anisotropy have been performed a con- 
sensus has not yet been reached. 

Several experiments indicate that the gap is finite at 
all points on the Fermi surface, however, others have 
suggested the presence of pronounced gap anisotropy or 
nodes. The presence of a fully gapped state is supported 
by point-contact tunneling spectroscopy measurements 
of SmFeAsOa;Fj^ (Sm-1111) P, [3, 0], penetration depth 
measurements of Pr-1111 Q^Sm-llll [|, Nd-1111 Q, 
and Ba2;Kj^Fe2As2 (Ba-122) [7] and angle resolved pho- 
toemission measurements of Ba-122 [1, M 11|. Sev- 
eral of these measurements show evidence of two distinct 
gaps. On the other hand, some experiments suggest the 
presence of low energy excitations which could be indica- 
tive of nodes. These include nuclear magnetic resonance 
measurements of La-1111 . [l2l .lT3| and penetration depth 
measurements of Ba-122 |14l |. 

Theoretically, there is considerable debate about the 
nature of the gap anisotropy and symmetry. Band- 
structure calculations have indicated that both the 1111 
and 122 materials have multiple quasi-two-dimensional 
sheets of Fermi-surface, hole-like close to the zone cen- 
ter and electron-like at the zone corner. At least for 
the undoped parent compounds these sheets are close 
to a nesting instability which can drive the system to- 
wards antiferromagnetism [l5| . Several papers have ar- 
gued that spin-fluctuations mediate the electron pairing 
and Mazin et al. il5] argued this favors an s± paring 
state where both the electron and hole Fermi surfaces 
are fully gapped but with gap-functions which are tt out 
of phase. Others however, have come to different conclu- 
sions. Most recently, Graser et al. [l^j have calculated 
that the two lowest energy gap functions, one with s sym- 



metry and one with d symmetry both have nodes on one 
or more of the Fermi surface sheets and are nearly de- 
generate. 

The iron phosphide superconductor LaFePO, is 
isostructural with LaFeAsO (La-1111). Its relatively low 
Tc ~ 6 K has been linked to the fact that the Fe-P 
bond angles depart substantially from those of a regu- 
lar tetrahedron Unlike the corresponding As-based 
compound, nominally undoped LaFePO is non-magnetic 
and superconducting [l^. The electronic structure of 
LaFePO has recently been explored in some detail by 
de Haas- van Alphen (dHvA) measurements [2^, which 
confirm the band-structure predictions of almost nested 
electron and hole pockets, a feature expected to be com- 
mon to all Fe-pnictides not subject to antiferromagnetic 
Fermi surface reconstruction. The detailed knowledge 
of the electronic structure and availability of supercon- 
ducting samples with very low levels of disorder makes 
LaFePO an ideal material in which to study order param- 
eter symmetry. In this letter, we report measurements of 
the London penetration depth of LaFePO down to very 
low temperature T < 0.02 Tc. Our data show strong ev- 
idence for the presence of line-nodes in this compound. 

Our samples were grown via a flux method which pro- 
duces plate-like single crystals with typical dimensions 
0.15 mm in the basal plane and <0.03mm in the inter- 
layer direction [l^ . The samples measured in this study 
are from the same batch as those used for dHvA mea- 
surements where the mean free path was estimated to be 
> 1000 A for the electron sheets and ~ 500 A for the hole 
sheets 20]. Heat capacity, measured using a modulated 
temperature technique, on the same single crystals used 
for the penetration depth measurements confirmed the 
bulk nature of the superconductivity. The data shown in 
Fig. [1] shows a sharp jump in C at Tc with the midpoint 
at ~5.6K with width -0.25 K. 

Measurements of the temperature dependence of the 
London penetration depth were performed with a high 
resolution susceptometer based on a self-resonant tunnel 
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FIG. 1: (a) Heat capacity anomaly in a LaFePO single crys- 
tal. AC is the difference between the heat capacity data in 
zero field and for B = 0.47 T (where Tc <2 K). (b) RF sus- 
ceptibility of the same sample close to Tc. 



diode circuit which was mounted in a dilution refrigera- 
tor. The circuit operates at ~ 14 MHz with an extremely 
small probe field {Hac < lOmOe) so that the sample is 
always in the Meissner state [2lj. Changes in the res- 
onant frequency are directly proportional to changes in 
the magnetic penetration depth as the temperature of 
the sample is varied. The calibration factor is determined 
from the geometry of the sample, and the total perturba- 
tion to the resonant frequency due to the sample, found 
by withdrawing the sample from the coil at low tempera- 
ture [2^ . The sample is mounted on a sapphire rod, the 
other end of which is glued to a copper block on which 
the Ru02 thermometer is mounted. The sample and rod 
are placed inside a solenoid which forms part of the res- 
onant tank circuit. The position of the sample within 
the solenoid could be varied in-situ allowing us to vary 
the RF field seen by the sample and hence check for RF 
heating which in-principle could cause the thermometer 
and sample to be out of equilibrium. The sapphire sam- 
ple holder is of very high purity and has a very small 
paramagnetic background signal which was difficult to 
determine precisely and has not been subtracted. It cor- 
responds to at most 20 A at the lowest temperature. 
Susceptibility measurements in this apparatus are shown 
in Fig. [1] in the region near Tc- The mid point of the 
transition occurs at 5.6/,K consistent with the heat ca- 
pacity measurement on the same sample. Other samples 
measured had Tc values in the range 5.4-5.9 K. 

Fig. [2] shows the temperature dependence of the in- 
plane penetration depth, AAa(T) in three samples, all 
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FIG. 2: Temperature dependence of the in-plane penetration 
depth, AXa in three single crystals of LaFePO. Data for sam- 
ple #3 has been multiplied by 1.5. The curves are offset for 
clarity. Solid lines are power-law fits giving an exponent of 
1.2 ± 0.1. Curve (a) is an extrapolation of the two gap be- 
havior found in Sm-1111 [^. Curve (b) is the temperature 
dependence expected from weak coupling BCS behavior (i.e. 
A = 1.76 Tc). We have used A(0) ~ 2500A in both cases and 
multiplied both curves by 50 to make them visible on this 
scale. 



measured with the ac field perpendicular to the conduct- 
ing planes so only in-plane currents are induced. The 
temperature dependence is very similar in all three sam- 
ples and is approximately linear. A variable power law 
fit, i.e. AA(r) cx T" for T < IK gives n = 1.2 ±0.1. The 
close to linear temperature dependence of A is strongly 
indicative of nodes in the order parameter. 

The temperature dependence of A arises from the de- 
pletion of the screening superfluid by thermal excitations. 
The superconducting energy gap, A^, dictates the spec- 
trum of these excitations; if A^ is finite for all momentum 
wavevectors, an exponential temperature dependence of 
A(T) is observed. By contrast, in a superconductor with 
an unconventional pairing symmetry the order parameter 
passes through zero for some wavevectors, leading to the 
presence of low energy excitations down to low temper- 
atures and a power-law temperature dependence of the 
penetration depth. In a clean superconductor with line 
nodes on a quasi 2D fermi surface the angular averaged 
density of states has a linear energy dependence, leading 
to a linear behavior of A(r) at low temperature. 

The linear temperature dependence in LaFePO is in 
marked contrast with that observed in conventional su- 
perconductors, or those displaying anisotropic s-wave su- 
perconductivity, such as MgB2 or NbSe2 [23|, , or the 
As-based 1111 compounds described in the introduction. 
In an s-wave multiband superconductor the low temper- 
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FIG. 3: Low temperature normalized superfluid density for 
samples #1 and #3 using A(0) = 2500A. Solid lines are linear 
fits to the data. Inset shows the superfluid density over the 
full temperature range. The dashed line shows the behavior 
found for Sm-1111 



ature behavior is dominated by the smallest gap on the 
Fermi surface. We show in Fig. [2] the expected tempera- 
ture dependence from gaps of 1.76 Tc (the BCS weak cou- 
pling value) and 1.0 Tc (similar to that found in Sm-1111 
[5|). These show a clear exponential saturation which 
is not visible in our data. We find that if there is a 
small residual gap it at least 30 times smaller than the 
gap maximum. The presence of nodes is therefore much 
more likely. 

The absolute values of AA(r) are somewhat variable 
between samples. As shown in Fig.[2l dX/dT ~ 350 A/Hz 
in samples #1 and #2 and smaller value of 220 A/Hz 
in sample #3. This discrepancy in slope is likely to be 
due to uncertainties in the calibration factors relating 
the measured frequency shifts to AA(T) arising from sur- 
face roughness. LaFePO has a mica-like morphology and 
samples #1 and ^2 had visibly more flaky edges than 
sample#3, which had mirror like surface on all faces (a 
similar effect was observed in NbSe2 Alternatively, 
it could result from intrinsic differences in the superfluid 
density as a function of doping i.e., slight changes in the 
oxygen stoichiometry (the samples had slightly different 

Te). 

In the usual weak-coupling theory for superconductors 
with line nodes, the superfluid density ps — [A(0)/A(T)]^ 
varies linearly over a much wider range of T than A(T) 
(corrections to the former are OT^ whereas to the later 
they are of OT^ ) . Muon smn relaxation measurements of 
polycrystalline samples [25| report a depolarisation rate 
of 1.2^m~^, which corresponds to a penetration depth of 
A(0) ~ 2400 A (we note that these authors indicate some 
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FIG. 4: Temperature dependence of the penetration depth, 
AAniix, measured with H \\ ah. Data for the same crystal with 
II c are also shown for comparison. 



uncertainty in this result due to a comparable anomalous 
contribution to the depolarisation rate in those samples). 
Measurements in the related materials La-1111 [26| and 
Sm-1111 [27| span the range 2000-3500 A. Choosing a 
representative value of A(0)=2500 A we calculate the su- 
perfluid density as shown in Fig. [31 The flgure shows that 
despite the uncertainties in the ratio AA(r)/A(0) the su- 
perfluid density varies almost perfectly linearly with T 
from below ~ 2K to around 100 mK. This key result is 
robust to the assumed value of A(0), so there is strong 
evidence for nodal excitations in the data regardless of 
this choice. In a nodal superconductor the low temper- 
ature behavior is dominated by the excitations near the 
gap nodes. For a superconductor with line- nodes A(r) 
is linear, with a slope determined by the rate at which 
the gap grows away from the nodes. In the simple d- 
wave case with Afe = Aq cos(26'), AA = In 2A(0)/cBr/Ao. 
For the weak coupling value Aq ~ 2.14 fc^Tc this pre- 
dicts a value of around 100-200 A/K (for A(0) in the 
range quoted above), somewhat smaller than the slope 
observed experimentally in sample #1 and #2, but close 
to the value seen in sample #3. As this calculation relies 
directly on A(0), better measurements of this parameter 
are required to make any precise comparisons to this or 
other candidate gap functions. 

Differences in the temperature dependence of Ac and 
Aa can sometimes reveal information about gap symme- 
try. For example, in multi-band systems the c-axis re- 
sponse can also be different due to the presence of sheets 

For 



with differing anisotropy and different gaps [24 , 
sample ^1 measurements were also performed with the 
field parallel to the conducting planes. This induces a 
combination of in-plane and inter-plane currents and thus 
probes a mixture of AAq and A Ac. Sample #1 has a ra- 
tio of thickness to width of about t/w ~ 7 which mixes 
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in a substantial component of Ac- We denote this mix- 
ture AAniix- This data is shown in Fig. |4] along with 
data from H||c for comparison. The changes in A,„ix 
are larger, with dAmix /dT = 1300 A/Hz. Although the 
present crystals are not of sufficiently regular geometry 
to isolate the changes in Ac precisely, the data are con- 
sistent with Ac(T') also having a linear temperature de- 
pendence but with a much larger slope due to substan- 
tial electronic anisotropy. Measurements of the upper 
critical field at low temperature indicate an anisotropy 
^^2/^02 — 11 consistent with the quasi-2D character 
of the Fermi surface [i^]. When line nodes are present 
on a simply warped quasi-2D Fermi surface, as found 
in LaFePO, Ac(r)/Ac(0) is expected to be similar to 
Aa{T)/Aa(0). This is observed in the quasi-2D organic 
superconductors, where the in-plane and inter-plane pen- 
etration depth have a similar temperature dependence 

The importance of the observation of a linear T de- 
pendence of the superfluid density over a wide range of 
temperature is that it is very difficult to produce this 
from extrinsic effects. The identification by Hardy et al. 
(30| of a similar predominately linear behavior of A(r) in 
YBa2Cu307 was instrumental in identifying the dx2_y2 
state in the cuprates. Impurities will weaken any intrin- 
sic linear temperature dependence, with A(T) varying like 
in the dirty limit (3ll |. However, power-laws with ex- 
ponent close to two can result frorn, for example, a fully 
gapped s± state with impurities [33] ■ The impurities 
produce a non-zero zero energy density of states simi- 
lar to that found in the d-wave case, and this has been 
shown to account well for the NMR data 3^ S^l- Hence, 
although power-laws (with n 2) point towards an un- 
conventional gap symmetry they can be produced either 
with or without intrinsic nodes. However, the linear de- 
pendence observed here strongly points towards intrinsic 
nodes. 

The observation of an exponential behavior of A(r) is 
also a robust feature which is not easily mimicked by ex- 
trinsic effects, so the different behavior found here for 
LaFePO compared to the As based 1111 compounds de- 
tailed above, probably points towards these compounds 
having intrinsically different gap symmetry / anisotropy. 
Just such a possibility was recently emphasized in the 
theoretical work of Graser et al. [ig] • 

After completion of this work we learnt of local scan- 
ning SQUID susceptibility measurements in LaFePO 
crystals which also show a linear temperature dependence 

From precise measurements of the penetration depth 
in clean crystals we have shown strong evidence to sug- 
gest that the ferro pnictide superconductor LaFePO has 
nodes in its superconducting gap function. Combined 
with the detailed information about the full three dimen- 
sional electronic structure of this material derived from 
dHvA measurements, this should help to decide between 



competing theories for the origin of superconductivity in 
these materials. 
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